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VAPOR INGESTION IN A CYLINDRICAL TAJIK 
WITH A 

CONCAVE ELLIPTICAL BOTTOM 
by Andrew Klavina 


DUCT ION 


The objective of thle study is to estimate the propellant residual in the 
Centaur vehicle. The propellant considered is liquid hydrogen; flow rates and 
acceleration levels of Interest are summarized in Table I, and Figure 1 shows 
details of the tank geometry. 

The residual estimates are based on the approximate analysis of Ref. 1, as 
extended to low - g conditions in Ref. 2. Due to the level of approximation 
involved in this analysis, the numerical results obtained must be interpreted 
with caution, and experimental support would be desirable. 

In this study, the flow rate and acceleration level are represented by the 
Weber number W and Bond number B. The group W/(1+B) is assumed to correlate 
vapor ingestion data at all W and BJ using the Bernoulli equation for a surface 
streamline, this group is expressed in terms of the liquid depth, the Bond 
number, and two functions that involve the fluid velocity field and surface shape. 
Plausibility arguments are used to approximate these functions in a manner 
appropriate to the Centaur geometry. These approximations are then used to 
generate numerical values of W/(1+B) as a function of the bulk liquid depth and 
Bond number, and thus arrive at quantitative estimates of the propellant residual 
at the conditions of interest. 
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ANALYSIS 

The essential features of the vapor ingestion analysis presented in Ref. 1 
and modified to low-g conditions in Ref. 2 consist of (1 ) writing the Bernoulli 
equation for a surface streamline between a point immediately above the drain 
and a point far from the drain, (2) expressing the fluid ‘velocity at the near 
point in terms of the flow rate, and the pressure drop along the streamline due 
to surface tension in terms of the surface curvature, and (3) applying a criterion 
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that defines incipient pull-through to relate the liquid depth at pull-through 
to flow rate and acceleration level. In this analysis, several approximate, 
limiting expressions are derived for the relations required in step 2. 


Referring to Figure 2, the Bernoulli equation between a near point n Q n and 
far point M *e H is 


p o + pgh o + pu o = p «» + * 


(D 


where p is the liquid pressure Just below the surface, Pis the liquid density, 
g the acceleration, directed as shown, h is the fluid surface height above a 
reference level, and u in the fluid velocity immediately above the drain. The 
velocity at «o is assumed to be zero. Definition of a control area A serves to 
relate to the volumetric flow rate Q, 


u o k *fh~ • 


( 2 ) 


and the pressure difference p - p Is related to the surface curvature J, 

go 0 

P« - P Q = - J 0 ), (3) 

where <r is the surface tension. Substituting (2) and (3) into (1) and introd- 

2 2 2 3 

ducing the Bond number B = pga A and the Weber number W = p Q /tt <xa , where i 
is the tank radius, the Bernoulli equation becomes 

w TlA = K - h - a 2 (J oe - J )/B. (4) 

B 2A 2 

In Ref. 3, vapor ingestion data are correlated in terms of the Froude number 
2 5 

F = Q /ga ? at high g and the Weber number at zero g; the analysis of Ref. 2 as 
well as current work at LMSC suggest that vapor ingestion data at all acceleration 
levels may be correlated in terms of the group W/(1+B), which has the appropriate 
limiting values at g— o, «. Accordingly, (4) will be rewritten in terms of this 
parameter group. 


It is the quantities A and - J Q that depend on the tank geometry and will 
be approximate In the present study; for convenience, the # surface curvature term 
is replaced by the dimensionless function f defined as 


h . “ h n 

g o 

a 2 (J te - J ) 
• o 


f 


(5) 



Tha Bernoulli equation then becomes 


1 +BL. (b 

1+B |n « 


this relation depends on the tank geometry only through the control area A and 
the surface curvature function f. 



At incipient pull-through, h Q decreases much more rapidly than h w ; this 
criterion ie introduced into the Bernoulli equation by differentiating with 
respect to h Q and setting d^/dh o = 0 at h Q = h QC , h w = l^s 

h - h = 

•oc OC L 

This equation relates h and h Q at incipient pull-through, denoted by the sub- 
script c. Setting h„ = h and h = h in (6) and using (7) to eliminate one 
of these variables yields the relationship between W/(1+B) and residual depth, 
with Bond number as a parameter. 


2 _§A _1 1 $£_ T 1 (7) 

A 3 h " 1+Bf f 3h 

0 o J c 


The results obtained so far are applicable to any tank geometry, subject 
only to the assumptions of irrotational, inviscid, quasi-steady flow, that the 
fluid fills the drain end of the tank, and that vortexing is not present. The 
analysis is specialized to a particular tank geometry through the functions A 
and f. 

The control area function A 


The choice of this function depends on the quantity of propellant in the 
tank and on the tank geometry. Several choices come to mind for various fluid 
levels, as detailed in Figure 3. 

If the fluid surface is far removed from the drain, the velocity potential 
is well approximated by that of a point sink, as discussed Ref. 2, and the 
appropriate control area is a spherical surface, centered near the drain. This 
situation is depicted in Figure 3a. In view of the approximate nature of the 
analysis, the complexity of an exact expression for this 'area is not justified. 
Rather, the approximations indicated in the figure are used. Taking the sphere 


3 



IHBC/D 30681*5 
TXA 2733 


radius to be the dip height h Q , with b the slml -minor axis of the ellipsoid that 
form the tank bottom, A(h 0 ) In this case Is given by 


where 


d/h 


A(h Q ) * th o I 1 - 1 /6(h o /a) - 2/3 J\ - (d/h 0 ) ] 
/l + (h o /b) ‘ <: (a i: /b ii - 1) - 1 


(h/a) (s 2 /b 2 - D 


( 8 ) 

(9) 


This result corresponds to the Halt of h Q0 much larger than drain size discussed 
in Ref. 2. 


If the fluid level is low enough, the velocity field will be essentially 
circumferential, and the control area Is planar, as shown in Figure 3b. 


The corresponding control area is 


A(h ) 
o 



(h /b) >/ 1 - (h /b)* + sin’ 1 (h /b) 
0 0 o 


(io) 


The surface curvature function f 

The expression to be used for this function depends strongly on the acceler- 
ation level, as observed in Refa. 4 and 5. At large Bond numbers, one expects an 
essentially flat liquid surface, while at small Bond number and a small contact 
angle, the surface is highly curved. As in the computation of the control areas, 
simple plausible shapes are assumed, as illustrated in Figure 4. 

Cases (a) and (b) refer to the situation h Q /b>1, where the control area is 
given by Equations (8) and (9). With B— 0, the liquid surface is highly curved, 
and an order of magnitude estimate gives 

f = 1/2, h /b >1, B« 1. (11) 

o 

With B-»«>, the Interface is essentially flat, and 

f = 1/4, h Q /b>1, B » 1 , (12) 
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Casts (o) and (d) rafar to tbs situation h Q /b < 1 , eorraspondlng to the 
control area of Equation (10). In the B— 0 liaiti the surface of tha bulk 
liquid is atlll curved as in oaaa (a) but in tha dip region the surface is 
visualised as a section of torus, with radius R 0 «R„ . 


The expression for f in this oasa is 


f(h 0 ) = 1/2 


1 - x/l - (h/b)^ 


, h o /b<1, B« 1 . (13) 


In the B— • liait of oase (d), tha liquid surfaoe is a flat ring, depressed fay 
the amount - h in the vleinity of tha drain, and f is given fay 


f = v 2 /8, h 0 / b<1 » B>>1 * 


(U) 


When these expressions are substituted into Equations (6) and (7), the value 
of W/(1+B) above vhioh pull-through will occur can be computed as a function of 
liquid depth. Such calculations are discussed below. It is emphasized that the 
numerical values are a coarse approximation at best, and empirical support is 
desirable, particularly in the vicinity of h^/b «1, where it is not clear 
whether the deep - or shallow - tank analysis is appropriate, and at low B, where 
a relationship between the fluid height at pull-through and residual volume must 
be assumed. 


fig 

Figure 5 is a summary of the results of calculations carried out using the 
various control area and surface ourvature functions assumed for limiting cases* 
Each curve shows the fluid height at which vapor ingestion may be expected to 
occur for a given value of W/(1+B). Alternately, the curves give the maximum 
value of W/(1+B) allowable at a given fluid level if vapor ingestion is to be 
avoided. The volumes indicated in Figure 5 at several pull-through depths assume 
that the liquid surface is flat, as in the B-»«e limit. 

e 

Referring to Table 1, an estimate oan be made of whether or not pull-through 

will occur at the flow rates and acceleration levels of Interest. During the 

v 

* The computet lam isqns n ce is swMarised in fisblm II. 
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start-up transient, with no settling thrust, W/(1+B) peaks at about 9.3; 

extrapolation yields a pull-through height h *>3.3, a value well above that 

corresponding to an initial propellant mass of 1467 lb . If a settling thrust 

of 24 11^ is used, the peak W/(1+B) is reduced to 6.3 x 10 and the pull-through 

height is h /b *1.1, The corresponding number for a settling thrust of 100 lb, 

•c * 

is h #c /b m .83. These values are obtained with the control area assumed spherical; 
the pull-through height for the 100 lfc^ thrust case assuming a planar control 
area- is h^/b *1,1. This difference represents an uncertainty of a factor of 
two in the residual volume calculation in the neighborhood of b^/b * 1 ; the 
actual case is expected to be somewhere between the two. From a design viewpoint, 
however, there is no problem if the initial level is at l\^/b * 1.56 corresponding 
to a 27.4$ propellant loading. 

An additior ' observation concerns operation of the vehicle at the full thrust 
of 30,000 lbf. As propellant is expended, the vehicle mass decreases and the 
acceleration increases, decreasing W/(1+B). The propellant level also decreases, 
however, lowering the maximum value of W/(1+B), and introducing the possibility 
of vapor ingestlou at full thrust. Figure 6 shows the actual and maximum values 
of W/(1+B) as a function of burn time for the full thrust case assuming an initial 
LH 2 volume of 350 ft , a constant LH 2 flow rate of 1241 gal/min, and several 
oxidiser/fuel ratios. On the basis of these calculations, one does not expect to 
encounter vapor ingestion before the bulk liquid level drops below the top of the 
sump region. 
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3 2 

Propellant t liquid hydrogen, a/p = 26.6 any sec 

Initial vehicle weight W q = 14,024 lb m 
Initial propellant loading - 27.4# » 350 ft 3 

Acceleration levola, baaed on W : 

o 

Thruat T ■ 24 lb f B = 1.47 x 10 3 Settling 

T = 100 lb f B = 6.13 x 10 3 Settling 

T = 30,000 lb f B = 1.84 x 10 6 Main Engine 

Flow rate a: 
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(l) Pick load nuaber and cboloe of A ud f. 

(U) 

(1U) Bvaluate A(M ) and f(l ) 

o o 

(lv) lvaluate the derivatives 

■§a * jaa [* <*„♦ .01) • a<« 0 )] , 

o 

11 „ ua [* <I„ ♦ .01) - «*„)] • 

O 

(v) Substitute the results of (ill) aod (iv) into lq. (7) 
to obtain (b /h_) _ 

•0 O C 

(vl) Substitute the results of (iii) sad (v) into lq. (6) 
to obtain W/( l+l) 
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pig, 3 M»p«oyi^Arrio>Ji to c.oJt£oL pukJ ctioU. 


(a) K 0 /t > 1 . 
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(L) t, ,/t>3 j *»< 
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